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The next-to-minimal supersymmetric model with a hght doublet-like CP-odd Higgs boson and 
small tan/3 can satisfy all experimental limits on Higgs bosons even with light superpartners. In 
these scenarios, the two lightest CP-even Higgs bosons, hi and h2, and the charged Higgs boson, /i"*", 
can all be light enough to be produced at LEP and yet have decays that have not been looked for or 
are poorly constrained by existing collider experiments. The channel hi aiUi with ai t'^t~ 
or 2j is still awaiting LEP constraints for m^j > 86 GeV or 82 GeV, respectively. LEP data may 
also contain e"'"e~ — > /i2ai events where /i2 Za\ is the dominant decay, a channel that was 
never examined. Decays of the charged Higgs bosons are often dominated by — > Ty*'*'ai with 
a\ —> gg,cc,T'^T~ . This is a channel that has so far been ignored in the search for t — > h'^b decays 
at the Tevatron. A specialized analysis might reveal a signal. The light ai might be within the reach 
of B factories via T —> 701 decays. We study typical mass ranges and branching ratios of Higgs 
bosons in this scenario and compare these scenarios where the ai has a large doublet component to 
the more general scenarios with arbitrary singlet component for the ai. 
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I. INTRODUCTION 

Discovery of Higgs bosons and exploration of their 
properties is the key to understanding electroweak sym- 
metry breaking and a major step in uncovering the ulti- 
mate theory of particle physics. The Higgs boson is the 
last missing piece of the standard model (SM). In the- 
ories beyond the SM the Higgs sector is typically more 
complicated, 6.17. in the minimal supersymmetric model 
(MSSM) there arc two Higgs doublets which lead to five 
Higgs bosons in the spectrum: light and heavy CP-cven 
Higgses, h and H, a CP-odd Higgs, A, and a pair of 
charged Higgs bosons, H^. In the next-to-minimal su- 
persymmetric model (NMSSM) which contains an addi- 
tional singlet superfield with complex component scalar 
field there are three CP-even Higgses, /ii.2,3, two CP-odd 
Higgses, ai.2 and a pair of charged Higgs bosons; and 
there are many simple models with an even more com- 
plicated Higgs sector. 

Since searches for Higgs bosons rely on detection of 
their decay products, it is crucial to understand the way 
the Higgs bosons decay. Although it is usually the case 
that one of the Higgs bosons has couplings to the W, Z 
bosons and to fermions that are close to those of the SM 
Higgs, it is not necessarily true that such a Higgs decays 
in the way the SM Higgs does [l[. A significant model 
dependence of decay modes applies to other Higgses as 
well. 

It has been recently argued that supersymmetric mod- 
els in the region of parameter space for which tan/3 is 
small and, in addition, there is a light doublet-like CP- 
odd Higgs predict that all the Higgses resulting from the 
two Higgs doublets {h, H, A and H^) could have been 
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produced already at LEP or the Tevatron, but would 
have escaped detection because the decay modes have 
either not been searched for (or the searches have been 
incomplete) or are ones to which the experiments are 
not sensitive. Although this scenario is ruled out in the 
MSSM, it is only marginally disallowed for niA < 2mi, 
and tan /? ^2.5 and thus can possibly be viable in sim- 
ple extensions of the MSSM The reason is that for 
niA <C mw and tan/? 1 the light CP-even Higgs bo- 
son becomes SM-like, and although it is masslcss at the 
tree level, its mass will receive a contribution from su- 
perpartners and the tree level relation between the light 
CP-even and CP-odd Higgses, to^ < rriA, is dramati- 
cally changed by SUSY corrections. Even for modest 
superpartner masses the light CP-cven Higgs boson will 
be heavier than 2mA- In particular, for superpartner 
masses between 300 GeV and 1 TeV and tan/3 ~ 1, one 
finds mil ~ 40 — 60 GeV and the h AA decay mode is 
open and generically dominant. 

Since the h has SM-like WW, ZZ couphngs, e+e^ 
hA is highly suppressed and the limits from the Z width 
measurements can be easily satisfied even for m.h + mA < 
mz- On the other hand, the e+e^ s- Zh cross sec- 
tion would be maximal. However, for tan/3 ~ 1 and 
mA < 2mb the decay width of the A is shared between 
T+r^, cc and gg final states and thus the (dominant) 
h AA decays are spread over many different final 
states: 4r, 2t2(7, Ag, 4c, 2g2c, 2r2c and 65, the latter 
being greatly suppressed relative to the SM expectation 
due to the presence of the h — > AA decays. As a con- 
sequence, the LEP limits in each channel separately are 
very substantially weakened. Of course, the decay mode 
independent limit requires a Higgs with SM-like ZZ cou- 
pling to be above 82 GeV It is this fact that rules 
out this scenario in the MSSM, since there mh cannot be 
pushed above 82 GeV by radiative corrections. 

The rest of the Higgs spectrum is basically not con- 
strained at all in this scenario. The heavy CP-even and 
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the CP-odd Higgscs could have been produced at LEP in 
e+e" — > HA but they would have escaped detection be- 
cause H dominantly decays to ZA - a mode that has not 
been searched for. Additional constraints are discussed 
in detail in Ref. 0. The charged Higgs is also very lit- 
tle constrained and up to ^ 40% of top quarks produced 
at the Tevatron could have decayed into charged Higgs 
and the b quark since the dominant decay mode for the 
charged Higgs W^^*''A with A — + cc, gg or t~^t~ 

would not have been separated from the the generic top 
sample ^ . In addition, pair production of a charged Higgs 
boson with the properties emerging in this scenario and 
mass close to the mass of the W boson could explain 
the 2.8(7 deviation from lepton universality in W decays 
measured at LEP [1] as discussed in 0]. 

The mass of the light CP-even Higgs is the only prob- 
lematic part in this scenario. There are however var- 
ious ways to increase the mass of the SM-like Higgs 
boson in extensions of the MSSM. A simple possibility 
is to consider singlet extensions of the MSSM contain- 
ing a XSHuHd term in the superpotential. It is known 
that this term itself contributes X'^v'^ siv? 2(3, where v = 
174 GeV, to the mass squared of the CP-even Higgs Q 
and thus can easily push the Higgs mass above the decay- 
mode independent limit of 82 GeV. Note that this con- 
tribution is maximized for tan /3 ~ 1 . In addition, it 
need not be the case that the light CP-even Higgs has 
full strength ZZ coupling, in which case the model- 
independent limit on m/j is reduced, while at the same 
time the H which carries the rest of the ZZ coupling can 
have mass above the LEP kinematic reach and/or decay 
to modes for which the LEP limit of 114 GeV does not 
apply. Thus, it is not surprising that in the NMSSM it 
is possible to find scenarios in which the lightest CP-odd 
Higgs has mass below 2mi, and the two lightest CP-even 
Higgs bosons and the charged Higgs would all have been 
produced at LEP and yet escaped detection. 

In this paper we study NMSSM scenarios with a light 
CP-odd Higgs boson and small tan/3. We will in par- 
ticular examine the subset of these scenarios in which 
the light CP-odd Higgs boson is mainly doublet-like 
(ai-doublet-like scenarios) and will find that they have 
many features in common with the MSSM scenarios dis- 
cussed above, except that they are not ruled out by 
Higgs searches — they are phenomenologically viable 
even with very light superpartners. For the subset of 
the fli-doublet-likc scenarios in which the hi has nearly 
SM-like couplings, the hi can be as light as 82 GeV (the 
decay-mode independent limit) by virtue of dominant 
decays hi — > aioi — > 2t2c, 4t, 4c, etc.. There are also 
scenarios for which the hi has reduced coupling to ZZ , 
dzzhi/dzzhsM ~ ^ these latter cases, the CP-even 



Higgs boson can have nihi as low as ~ 55 GeV. All these 
scenarios are similar to the scenario with a light singlet- 
like CP-odd Higgs in the NMSSM 0, i, i, M EM in 
that it is the unexpected Higgs decays that allow one or 
more light Higgs bosons to have escaped LEP detection. 
The irnportant difference is that the scenario discussed in 
Refs. 0, H, [^, 113, El is the usual decoupled scenario as 
far as the two Higgs doublet part of the Higgs spectrum is 
concerned (the CP-odd Higgs, the heavy CP-even Higgs 
and the charged Higgs are heavy and approximately de- 
generate) and the light CP-odd Higgs is supplied by the 
additional singlet. In contrast, in the oi-doublet-like low- 
tan /3 scenarios, the CP-even and CP-odd Higgses coming 
from the additional singlet are typically heavy and do not 
drastically alter the two Higgs doublet part of the Higgs 
sector. The latter then looks like the Higgs sector of the 
MSSM with somewhat modified mass relations. 



II. LIGHT DOUBLET-LIKE ai IN THE NMSSM 

As already mentioned the scenario with a light doublet- 
like CP-odd Higgs and small tan /3 is phenomenologically 
viable in the simplest extension of the MSSM, the next- 
to-minimal supersymmctric model which adds only one 
singlet chiral superfield, S. The very attractive nature of 
the NMSSM extension of the MSSM on general grounds 
has been discussed for many years [T^ : in particular, it 
avoids the need for the fi parameter of the MSSM super- 
potential term ^HuHd- The NMSSM particle content 
differs from the MSSM by the addition of one CP-even 
and one CP-odd state in the neutral Higgs sector (assum- 
ing CP conservation), and one additional neutralino. We 
will follow the conventions of 15]. Apart from the usual 
quark and lepton Yukawa couplings, the scale invariant 
superpotential is 

A SHuHd + ^ ^3 (1) 

depending on two dimensionless couplings A, k beyond 
the MSSM. [Hatted (unhatted) capital letters denote su- 
perfields (scalar superfield components).] An effective /i 
term arises from the first term of Eq. ([T]) when the scalar 
component of S acquires a vacuum expectation value, 
s = {S), yielding 

/icff = As . (2) 

The trilinear soft terms associated with the superpoten- 
tial terms in Eq. ([1]) are 

XAxSH^Hd + ^A^S^ . (3) 



^ We thank Ricardo Euscbi (CDF) for a detailed discussion of the 

CDF and DO analyses. plings of the Higgses to Z and W through mixing, see e.g. 

^ In general, in singlet extensions it is possible to alter the cou- Refs. \Ti Il3f or to provide new Higgs decay modes [y. 
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The final input parameter is 

tan/3 = , (4) 

where /i„ = {H^}, hd = (Hd). The vevs hd and s, 
along with mz, can be viewed as determining the three 
SUSY breaking masses squared for Hd and S (de- 
noted 77ijj , iTijj^ and m|) through the three minimiza- 
tion equations of the scalar potential. Thus, as com- 
pared to the three independent parameters needed in the 
MSSM context (often chosen as /i, tan/3 and Ma), the 
Higgs sector of the NMSSM is described by the six pa- 
rameters 

A , K , , tan/3 , ^off • (5) 

(We employ a convention in which all parameters are 
evaluated at scale mz unless otherwise stated.) We will 
choose sign conventions for the fields such that A and 
tan/3 are positive, while k, A\, A^^ and ^off should be 
allowed to have either sign. In addition, values must be 
input for the gaugino masses (M\;2..-i) and for the soft 
terms related to the (third generation) squarks and slep- 
tons (toq, m|), m^^, m|;. At, At, and A^-) that con- 
tribute to the radiative corrections in the Higgs sector 
and to the Higgs decay widths. For small tan/3, the soft 
parameters which play the most prominent role are ttiq, 
roij and At- 

A complete survey of the parameter space is difficult. 
To present results in a manageable way, we fix /i and 
tan j3 together with all soft SUSY breaking masses and 
scan over trilinear and soft-trilinear couplings. We will 
plot results in various two-dimensional planes. The input 
parameters of Eq. ^ are scanned over the following re- 
gions with fixed steps: A e (0.001,0.6) using 60 steps of 
size 0.01; k G (—0.6, 0.6) using 120 steps of size 0.01, with 
some refined scans for k G (—0.06, 0.06) using 120 steps 
of size 0.001; Ax E (-600 GeV, 600 GeV) using 200 steps 
of size 6 GeV; and finally A^ G (-600 GeV, 600 GeV) us- 
ing 200 steps of size 6 GeV, with refined scans for this 
same range with 1000 steps of size 1.2 GeV. Varying the 
fixed soft SUSY breaking masses leads to smaller changes 
than does varying tan /3. Thus, we will consider only a 
few choices of soft SUSY breaking masses and will focus 
on the important changes that occur as tan/3 is changed. 

All scans are performed in the context of NMHDE- 
CAY. NMHDECAY checks a long list of experimental 
constraints, especially those coming from LEP data. It 
also checks various theoretical constraints on the model 
incorporated, such as requiring that the vacuum be a true 
vacuum. NMHDECAY also issues a warning if any of the 
couplings, A, K, ht or hi, become non perturbative {ht and 
/if, are the Yukawa couplings) after evolution to the GUT 
scale. We will consider scenarios in which these become 
non-perturbative as well as scenarios in which they re- 
main perturbative. Aside from this, all plotted points 
are consistent with all the NMHDECAY constraints. 

As stated in the introduction, we wish to focus on cases 
for which iriai < 10 GeV. Such scenarios have the most 



unusual features. In this mass region, it is important 
to incorporate the constraints arising from recently im- 
proved limits on B{T 701) with oi — > rr from CLEO- 
III as well as old CUSB-II limits HH on S(T ^ 701) 
where ai is only assumed to be visible. These basically 
place an upper limit on the bbai coupling defined by 

in the region TOqi < my. Further constraints on this 
coupling were obtained at LEP by looking for bbai pro- 
duction with ai — > rr and ai — > bb [2^ . [23| . The former 
channel is important in the < 10 GeV that we focus 
on. The upper limits on C^^;,^ using the above inputs 
arc given in Fig. 1 of [23|. At any given tan/3, a limit 
66 1 66 1 converts to a limit on | cos 6*^1 ViS- 
ingC„,66 = cos0Atan/3, i.e. |cos0a| < IC^^^sl"^""/ tan/3. 
The resulting values for \ cos9a\™'^^ appear in Fig. 3 of 
[23 |. It also turns out that Tevatron limits on pp 
ai M~''A'~ provide some constraints on C^^f,^ in the re- 
gion from 8 GeV < mg, < 9 GeV that are stronger than 
those from e+e" data [2^. These too are incorporated. 

A final addendum to NMHDECAY is to include off- 
shell decays involving an ai and a gauge boson in the 
final state. In particular, /i+ W*ai and /12 — * Z*ai 
virtual decays are of occasional importance in the small 
TOai region. 



A. Results for tan/3 = 2 

A convenient reference scenario is the case of tan /3 = 2 
with MsusY = 300 GeV and At ^ Ab ^ Ar ^ 
-300 GeV. The plots for this case are Figs.[I]-[ll In our 
plots, the blue -f 's are all points that satisfy the NMHDE- 
CAY constraints, while green diamonds are those which 
in addition have a light CP-odd Higgs which is doublet- 
like, cos^ 6a > 0.5. The red crosses single out those 
points for which mh^ < 65 GeV. Because of the lim- 
its on I cos^yil discussed above, values below about 
7.5 GeV arc disallowed for cos^ 9a > 0.5, as are many 
points with cos^ 9a < 0.5. This is illustrated in Fig. [T] 
The jagged shape of the boundary in the ma^ < 7.5 GeV 
region for the cos^ 6a < 0.5 points simply reflects the 
rather rapid variations in the limits from i3(T — > 701) 
decays. 

The first plot, Fig.[5]shows the m/^ masses that are ob- 
tained in our scan and the correlation with ma^ . Higgs 
with mhi < 114 GeV are not excluded by LEP data. 
The reason is apparent from Figs. O [3] and \5\ There, 
we plot Clihi) = g%zhj9lzhs.m^ B{hi aioi) and 
Cy[hi)B{hi bb) vs. m/^. We see that the light hi es- 
capes LEP constraints mainly because of large B{hi 
aiQi) (where oi — *■ t~^t~ yields a 4r final state that is 
weakly constrained by LEP data) although there are a 
significant number of points for which the ai is mainly 
singlet with small Cy(hi). The plot of Fig. [5] shows the 
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tan/S=2.0, ;ti=150 GeV Ms„sy=300 GeV, At=-300 GeV 
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FIG. 1: cos^ 6a is plotted vs. niai for the tan /3 = 2, msusY ~ 
300 GeV, A = -300 GeV scenario. 



FIG. 4: B{hi — > aifli) is plotted vs. rrih-^ for the tan/3 = 2, 
msuSY = 300 GeV, ^ = -300 GeV scenario. 
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FIG. 2: TTiai is plotted vs. for the tan/3 = 2, msuSY 

300 GeV, A = -300 GeV scenario. 
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FIG. 5: Cv{hi)B{hi bb) is plotted vs. for the tan/3 
2, msusY = 300 GeV, A = -300 GeV scenario. 



net rate for e~^e~ Zbb relative to the SM prediction. 
We observe that away from the 90 GeV to 105 GeV win- 
dow in ruhi = "^bbi ill which there is an excess of LEP 
events relative to background, this net rate must be quite 
small. In the 90 GeV to 105 GeV window, the best fit to 
the experimental data corresponds to an an excess of or- 
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FIG. 3: Cy{hi) is plotted vs. for the tan/3 — 2, msusY 
300 GeV, A = -300 GeV scenario. 



dcr 0.1 times the expected SM rate is allowed. However, 
in this window, an excess as large as 0.2 times the SM 
rate is still allowed at 90% CL, as reflected in the plot. 
Note that it is mainly the points with cos^ 6 a ^ 0.5 that 
best explain the observed 0.1 x SM excess in this region. 

The next interesting feature of these small rUai , small 
tan/3 scenarios is the very substantial probability that 
the will also be quite light. As shown in Fig. [51 this 
is particularly the case for parameters such that the oi 
is mainly doublet. For these ai-doublet-likc scenarios, 
we observe that there are cases for which m/^ is well 
below 100 GeV while m/j+ is of order 100 GeV, and 
in the vast majority of these ai-doublct-like scenarios 
< 170 GeV so that the would have been pro- 
duced in top decays. At the same time, as shown in 
Fig. [71 for the ai-doublet-like scenarios m/jj can also be 
of order 100 GeV, and in nearly all cases m/j^ < 200 GeV 
so that e+e^ Zh2 production events would be present 
in LEP data. In the case of ~ 100 GeV it is the re- 
duced Cy{h2) (Fig. [S]) coupled with large B{h2 — > aiOi) 
(Fig. [9]) that makes LEP sensitivity in the Zh2 Zbb 
channel small. Indeed, Fig. [TOl shows that the /12 con- 
tribution to the Zh2 — > Zbb channel can only be sig- 
nificant for ?Ti/i, ^ 125 GeV, well above the LEP kine- 
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FIG. 6: is plotted vs. m;,j for the tan/3 = 2, msuSY 

300 GeV, A = -300 GeV scenario. 



FIG. 8: Cl{h2) is plotted vs. mh,2 for the tan/3 = 2, msusY 
300 GeV, A = -300 GeV scenario. 
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FIG. 7: m;,2 is plotted vs. for the tan/3 = 2, msusY ~ 

300 GeV, yl = -300 GeV scenario. 



FIG. 9: B{h2 — > aifli) is plotted vs. for the tan/3 = 2, 

msusY = 300 GeV, A = -300 GeV scenario. 



matic reach. (However, as we shall see, this conclusion 
does not apply to all choices of tan (3 and soft-SUSY- 
breaking parameters.) For these same scenarios with 
large cos^ 9a, the ai and /i2 both have substantial doublet 
component, and the Z — > h2ai rate at LEP would also 
have been significant. For mji^ near 100 GeV, the /i2ai 
final states would have escaped LEP detection because 
of large i?(/i2 — > aiOi). For larger up near 200 GeV, 
B{h2 — > Zai) would have been large, see Fig. [TTl and 
LEP did not analyze their data in such a way as to be 
sensitive to /i2ai Zaiai final states, especially given 
that ai decays to either two taus or two jets. 

Turning to the charged Higgs boson, most LEP 
searches for the /i+ were based on the dominant e+e^ ^ 
h'^h~ production mechanism assuming that /i"*" —f t'^Vt 
and cs were the only two decay modes. However, 

Fig.[T2lshows that h+ W+ ai is dominant for the ai- 
doublet-like scenarios. Limits on mfi+ weaken as B{h~^ — > 
t~^i^t) declines — if B{h^ t'^Vt) ^ 1,0.5,0 the lim- 
its are roughly m;j+ > 90,80,80 GeV, respectively j26j . 
DELPHI considered the possibility of — > W'^'^*'^a as- 
suming a — > 66 is dominant [2^ . However, their limits 
on mh+ do not apply to the case of ma < 2mi, of interest 
here. 



Overall, we have the remarkable result that for the 
chosen tan/3 = 2 and tosusy = 300 GeV, A = -300 GeV 
parameters there are a large number of model points (the 
ai-doublet-likc points) for which the hi and have mass 
at or below 100 GeV and the /i2 has mass in the range 
100 — 190 GeV. All would have been copiously produced 
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FIG. 10: Cv{h2)B{h2 bb) is plotted vs. mn^ for the 
tan/3 = 2, msusY = 300 GeV, A = -300 GeV scenario. 
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FIG. 11: B{h2 —> Zai) is plotted vs. m^^ for the tan/3 — 2, 
msusY = 300 GeV, A = -300 GeV scenario. 



FIG. 12: B{h+ Vy+'*'ai) is plotted vs. 771^+ for the 
tan/3 = 2, msusY = 300 GeV, A = -300 GeV scenario. 
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at LEP, and yet all would have escaped LEP detection. 

The primary sensitivity of the Tevatron to the ai 
doublet-like scenarios with a light is through searches 
for ti production with one t decaying via t — > h~^b 
[29l . [sot ■ The recent preliminary Tevatron analyses 2^, 
set limits on the B{t h^b)B{h^ t^Vt) as a func- 
tion of mft+ . Fitting simultaneously the branching ratio 
product and (j{pp —>■ tt), the limit for mfi+ = 80 GcV 
is B{t h+b)B{h+ r+Vr) < 0.12. In the present 
scenario, all these searches would have suppressed sen- 
sitivity for the cases where m/j+ ~ 100 GeV. The rea- 
son is that /i+ W~^^*'>ai always has branching ratio 
> 0.5, and the oi decays primarily to r+r" for niai > 
2TO7. and to various lighter final states if rriai < 2771,-. 
B{h+ VF+(*)ai) is shown in Fig. [TH From Fig. [HI 
we see that B{h+ r+Vr) ~ 1 - B{h+ W^+(*)ai) for 
TO/J+ < rat + rnb- 

Of course, B{t — *■ hr^h) is tan/3 dependent. For 
tan/3 = 1.2, 1.7, 2 it is of order 0.3, 0.173, 0.126 for ■mh+_ ~ 
90 GeV. For mf^+ > mt + one finds that tb is 

the dominant decay. At tan/3 = 2, B{h+ VK+(*)ai) 
can be of order 0.5 out to relatively large m/j+. LHC 
search strategies sensitive to all these unusual scenarios 
need to be developed. Some discussion of the possibilities 
appears in Ref. [3l|. 

Given the fact that the ai appears in so many decays, 
it is useful to review its branching ratios. In NMHDE- 
CAY, these are computed using partonic final states and 
masses. This implies a few inaccuracies. In particu- 
lar, ai — > ss is non-zero even when < 2mK since 
NMHDECAY employs = 0.19 GeV. The ai branch- 
ing ratios appear in Figs. fMlfTSlfTBl andflTl As expected, 
if < 2'mb but well above 2?7i,t- (as is the case for all 
fli-doublet-like scenarios), oi — !• T+r~ is the dominant 
decay, with the remainder being in the oi — > gg and 
Oi cc modes (in that order). For cases where the ai is 
approaching 2™,-, B{ai t+t~) declines, but is always 
bigger than 0.5 if > 2mr with the residual mainly 
taken up by B{ai cc). For the few < 2mc points, 
ai ^ ss is dominant. 

As noted earlier, an important constraint on scenarios 
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FIG. 13: B{h+ -> T+Ur) is plotted vs. for the tan /3 = 2, 

?TisusY ~ 300 GeV, A = —300 GeV scenario. 



with a light ai is the branching ratio B(T 701 ) w hich 
is strongly constrained by data from CLEO-III [20| for 
TOqj > 2tot-. In particular, depending upon the precise 
value of rria^ in the range between 2mr and 7.5 GeV, the 
95% CL upper limit on B{T — > 701) is between 6 x 10~^ 
and 1.2 x 10"'"'. In Fig. [18] we plot B{T 701) as a 
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FIG. 14: B{ai — > r"'"r~) is plotted vs. niai for the tan/3 — 2, 
msusY = 300 GeV, A = -300 GeV scenario. 
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FIG. 15: B{ai —r cc) is plotted vs. niai for the tan /? = 2, 
msusY = 300 GeV, A = -300 GeV scenario. 



FIG. 18: B{T 701) is plotted vs. niai for the tan/3 = 2, 
msusY = 300 GeV, A = -300 GeV scenario. 
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FIG. 16: B{ai — > ss) is plotted vs. for the tan/3 = 2, 

msusY = 300 GeV, A = -300 GeV scenario. 



FIG. 19: j/max is plotted vs. A for tlie tan/3 — 2, msusY 
300 GeV, A = -300 GeV scenario. 



function of rUai after imposing B{T — > 701) constraints. 
We see that points with cos^ 6a > 0.5 (the green dia- 
monds) are only allowed at relatively large rUai ^^'^ that 
even some points with cos^ 9a < 0.5 have been elimi- 
nated in the < nia^ < 7.5 GeV range. Thus, it 
is the T 7ai — > 7r+T^ decay limits that rule out 
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FIG. 17: B{ai — » gg) is plotted vs. iriai for the tan [3 = 2, 
msusY = 300 GeV, A = -300 GeV scenario. 



ai-doublet-like scenarios with nia^ < 7.5 GeV. The un- 
derlying reason for the oi-doublet-like points to be more 
strongly excluded by T decays is that, as discussed ear- 
lier, the aibb coupling is given by C^^^^ = cos 0^ tan/?, 
which is, of course, largest for large | cos 6*^11. 

Let us turn for a moment to a discussion of whether 
or not these scenarios are perturbative after evolution to 
the GUT scale. The couplings of interest are A, k, ht 
and hh. At low tan/?, hb always remains perturbative 
but A, K and ht can become large. In Fig. [12] we plot the 



^ For small tan /3 the the top Yukawa coupling becomes non- 
pcrturbative close to the grand unification (GUT) scale. The 
exact value of tan /3 consistent with pcrturbativity all the way to 
the GUT scale depends on superpartner masses through SUSY 
threshold corrections to the top Yukawa coupling, and in the 
NMSSM it is about tan/3 > 1.6. However, adding extra vector- 
like complete SU(5) matter multiplets at the TeV scale, e.g. the 
parts of the sector that mediate SUSY breaking (messengers) or 
are present for no particular reason, does not affect the unifica- 
tion of gauge couplings while it slows down the running of the top 
Yukawa coupling [la . [ITII and even tan /3 ~ 1 can be consistent 
with perturbative unification of gauge couplings. 
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FIG. 20: 1/max is plotted vs. rimax for the tan/3 = 2, msusY ~ 
300 GeV, A = -300 GeV scenarfo. 

value of 

^ inax{X, K,ht,hb} , , 

?/max — [i ) 

at the GUT scale as a function of A for the various sce- 
narios in our tan/3 = 2 scan. A value of j/max = 0.5 
indicates that the evolution has gone non perturbative. 
In Fig. [201 we show which of the couplings is largest or 
has gone non-perturbative first using the code 4 = A, 
5 = K, 6 = ht and 7 = h},. We observe that it is most 
often K that has the largest coupling at the GUT scale, 
especially for the cos^ 9a > 0.5 scenarios. 

Another issue of interest is whether finetuning of the 
NMSSM parameters (in particular A\ and A^) is re- 
quired (either at scale mz or at the GUT scale) in order 
to obtain rUai < 10 GeV and scenarios that escape LEP 
and other limits. In we developed a measure G of 
this fine tuning. In Fig. [5T1 we plot G as a function of 
cos 9a ■ We see that small values of G arise for quite spe- 
cific values of cos 0a, namely — 0.6 < cos 6* a < — 0.4 and 
0.15 < cos9a ^ 0.22. Note that the ai-doublet-like sce- 
narios typically have moderately large G values — only if 
the fli is singlet like is it possible for there to be no need 
for tuning Ax and A^ in order to achieve < 2mh and 
large B{hi aiUi) (to escape LEP limits) simultane- 
ously. 

In our scans, we did not specifically exclude scenar- 
ios because of difficulties with precision electroweak con- 
straints (mainly the parameter T) or the anomalous mag- 
netic moment of the muon, a^. In fact, for all the points 
plotted, -0.002 < AT < 0.011, where AT is defined 
relative to a SM-like Higgs with mass 100 GeV, and 
-2.2 X 10-1° < da^ < -1.4 x lO'^" where Sa^, is the 
net contribution of the entire Higgs sector. Clearly, the 
size of AT is such that the Higgs sector of the NMSSM 
models being considered makes a very small contribution 
to T while Ja^ is also so small as to have little impact on 
the current discrepancy between the SM prediction for 
and the experimental observation, which difference is 
of order 30 x lO^^". 

Finally, we show in Figs. and the A, k, A\ and 



FIG. 21: G is plotted vs. cos^a for the tan/3 — 2, msuSY = 
300 GeV, A — —300 GeV scenario. The displayed points 
comprise only a small fraction of the total number of points 
appearing in previous figures. 
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FIG. 22: A is plotted vs. k for the tan /3 = 2, msusY = 
300 GeV, A = -300 GeV scenario. 

Af^ values which yield the points plotted in the preceding 
figures. The main observation is that the cos^ 9a > 0.5 
points require k > and Ak,,A\ < 0. Note also the small 
number of points with k close to zero. Many, but not all, 
of the very small rrih^ < 65 GeV scenarios arise from 
these points. 



B. Results for tan/3 = 1.7 

One can avoid non-perturbative couplings for a large 
number of allowed points for lower tan (3 if tosusy and At 
are somewhat larger. As an example, we present results 
for the case of tan/3 = 1.7, msuSY = 500 GeV and A ~ 
-1000 GeV in Figs, [^illi^ The point notation is as for 
tan/3 = 2, except that in this case there are points for 
which B(hi — > oioi) < 0.7. These points are indicated 
by the yellow squares. 

As in the previous case, significant restrictions are 
placed on | cos 0^1 due to limits on the C^^^^ coupling. 
Fig. [M] shows that once again these restrictions basically 
imply a limit on cos^ 9a that is significantly below 0.5 if 
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FIG. 23: Ax is plotted vs. Ai^ for the tan/3 — 2, msuSY ~ 
300 GeV, A = -300 GeV scenario. 



FIG. 26: Cv{hi) is plotted vs. for the tan/3 = 1.7, 

msusY = 500 GeV, A = -1000 GeV scenario. 
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FIG. 24: cos^ ^.4 is plotted vs. for the tan/3 

msusY = 500 GeV, A = -1000 GeV scenario. 



1.7, 



ma, < 7.5 GeV. 

We now repeat the same set of figures as in the tan f3 = 
2 case. Many of the same comments apply. Where ap- 
propriate we shall comment on differences. 

We note that we have found many more points with 
quite low nih, for this case as compared to the tan/3 = 2 
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FIG. 25: niai is plotted vs. for the tan/3 = 1.7, ttisusy 
500 GeV, A = -1000 GeV scenario. 



TABLE I: Selected points for which mh-, and corresponding 
he within the LEP excess region and the corresponding 
Clihi)B{hi bb) and Cv{h2)B{h2 bb) values. 





Cl{hj)B{hi - 


bb) 




Cl{h2)B{h2 ^ bb) 


93.1 


0.0684 




96.2 


0.1590 


90.7 


0.0560 




96.6 


0.1726 


90.2 


0.1171 




97.2 


0.1468 


88.3 


0.0557 




97.0 


0.1803 


87.8 


0.0974 




97.5 


0.1609 


90.7 


0.0560 




96.6 


0.1727 


92.7 


0.1748 




97.2 


0.1037 


90.9 


0.0599 




97.1 


0.1416 



scenario previously considered. Of particular interest is 
the fact that there are a significant number of model 
points for which m^^ is near 90 GeV and m^, is just 
below 100 GeV with Cl{hi)B{hi -> bb) ~ 0.1 -0.2 and 
Cy{h2)B(h2 bb) - 0.1 - 0.2. A particular subset of 
these can be identified in Figs. [28] and [33] as the yel- 
low squares with the above attributes. (However, there 
arc quite a few blue points that also satisfy these cri- 
teria.) The precise masses and CyB{h bb) values 
of the yellow-square points are tabulated in Table IIIBI 
These points appear in Fig. [M] as the yellow squares 
with cos^6'yi (more precisely cos^^i 0.1) and 
TOai ~ 4 GeV or 8 GeV. They appear in Figs. [JT] and 
l^as the yellow square points with k e [—0.046, —0.041], 
A - 0.14-0.15, - 6-7 GcV and Ax - 486-492 GeV. 
The cos 9 A, k and ^„ values indicate that these points 
are ones that are close to the Peccei-Quinn symmetry 
limit of the NMSSM. 

The reason that these points are of particular inter- 
est is that the two Higgs bosons combine to nicely ex- 
plain the LEP excess seen throughout the entire to^j^ £ 
[88 GeV, 100 GeV] mass region in the e+e~ Zbb chan- 
nel. The level of this excess corresponds to CyB{h 
bb) ^ 0.1 — 0.2 for any single h mass in this region. The 
masses of the two Higgs bosons are typically separated 
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FIG. 27: B{hi — + ffliai) is plotted vs. rrih^ for the tan /3 — 1.7, 
msusY = 500 GeV, A = -1000 GeV scenario. 
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FIG. 29: is plotted vs. m^.j for the tan/? — 1.7, msusY 
500 GeV, A = -1000 GeV scenario. 
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FIG. 28: Cv{hi)B{hi ^ bb) is plotted vs. mn^ for the 
tan/3 = 1.7, msusY = 500 GeV, A = -1000 GeV scenario. 
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FIG. 30: mh, is plotted vs. rrihi for the tan/3 = 1.7, msusY 
500 GeV, A'^ -1000 GeV scenario. 



by about 2ares where ares ~ 3 GeV was the LEP mass 
resolution. As a result, the combination of the two Hig- 
gses would give the broad excess observed. The manner 
in which Cyi3(ft. — + 66) ~ 0.1 — 0.2 is achieved is quite 
different for hi vs. ft,2- In the case of the hi, Cy{hi) (see 
Fig. is small and B{hi bb) is fairly large (because 
B{hi aiGi) is relatively small (see Fig. ET)) . In the case 
of the /i2, Cy(/i2) (see Fig. |3T)) is large and B{h2 66) 
is fairly small (because B{h2 ^ aiai) is relatively large 
(see Fig.|32]). 

These special points are also rather attractive in that 
they arc ones for which the couplings remain perturba- 
tive after evolution evolution to the GUT scale (see the 
yellow-square points with small A in Fig. [38|) . However, 
for these points the Af^-Ax fine-tuning measure G (see 
the yellow-square points in Fig. [40] with cos 9 a ~ 0.1) is 
somewhat large. 

As regards the a i -doublet-like (green diamond) points, 
we observe from Fig.[30]that the lower bound on m/i^ has 
been pushed to about 110 GeV vs. the ~ 100 GeV value 
obtained for tan/3 = 2. This means that Zh2 produc- 
tion at LEP would have been minimal or absent for such 
cases, but there would still have been a significant rate 
for /i2ai production. However, to repeat, LEP did not 



look for the relevant /i2ai ^ aiaiOi or /i2ai ZaiUi 
final states that would have been dominant (Figs. [5^ and 
[34l respectively). 
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FIG. 31: C^(fe2) is plotted vs. 171^2 for the tan (3 
msusY = 500 GeV, A = -1000 GeV scenario. 
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FIG. 33: Clr{h2)B{h2 hb) is plotted vs. m^^ for the 
tan/3 = 1.7, msusY = 500 GeV, A = -1000 GeV scenario. 



FIG. 36: B{h'^ — > t^i't) is plotted vs. m,j+ for the tan/3 
1.7, msusY = 500 GeV, A = -1000 GeV scenario. 



C. Results for tan/3 = 1.2 

We have also performed a scan for the case of tan [3 = 
1.2, TOsusY = 500 GeV and A = -1000 GeV. For the 
most part, results are very similar to those for tan/3 = 
1.7, TOSUSY = 500 GeV and A = -1000 GeV. One dif- 
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FIG. 34: B{h2 — > Zai) is plotted vs. for the tan/3 = 1.7, 
msusY = 500 GeV, A = -1000 GeV scenario. 



ference arises because the coupling of the oi to bb, pro- 
portional to cos 9a tan /3 is weaker for the lower tan f3 
value. This implies that the experimental upper lim- 
its on this coupling are less restrictive at a given value 
of cos9a- The result is that jcos^Ap > 0.5 is possi- 
ble for TOqj < 2to^ (i.e. below the values for which 
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FIG. 37: B{T yai) is plotted vs. rUai for the tan/3 — 1.7, 
msusY = 500 GeV, A = -1000 GeV scenario. 
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FIG. 38: j/max is plotted vs. A for the tan/3 = 1.7, msuSY ~ 
500 GeV, A = -1000 GeV scenario. 



tan/S=1.7, ^(=150 GeV Msusy=500 GeV, At=-1000 GeV 




FIG. 39; J/max is plotted vs. Timax for the tan/3 — 1.7, 
msusY = 500 GeV, A = -1000 GeV scenario. 

the aibb coupling is so strongly limited by CLEO-III re- 
sults). This is made apparent by comparing Fig. 1431 to 
Fig. [211 Another difference is that for all but a special 
class (to be described later) of the tan/3 =1.2 scenarios, 
one or more of the couplings, A, k, A\ or becomes 
non-perturbative in evolving to the GUT scale. 

In Fig. |44] we plot mh+ vs. nihi- We see that the 
lowest value of is about 90 GeV and arises for the 
ai-doublet-like scenarios. There arc a significant number 
of points with ruh^ < 65 GeV, all of which have a very 
singlet-like oi, as is most apparent from the red crosses 
in Fig. 221 Fig. [ISl shows vs. m^i . We see that 
nih2 values as low as 90 GeV are possible for singlet-like 
ai, whereas the lower limit on m/i^ for oi-doublct-like 
scenarios has risen to about 140 GeV as compared to the 
lower values found for tan/3 =1.7 and 2. 

For this tan/3 = 1.2 case, there are many points 
with TO/J2 in the 95 — 100 GeV interval and in the 
90 - 96 GeV interval with both C'^{hi)B{hi bb) and 
C^{h2)B{h2 -> bb) between 0.05 and 0.15 that would 
explain the broad excess in this region seen at LEP. 
As for tan/3 = 1.7, all these points have very small 
K and An and are therefore close to the PQ symmetry 
limit of the model. Most of these points are such that 
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FIG. 40: G is plotted vs. cos 6a for the tan/3 = 1.7, msuSY = 
500 GeV, A = -1000 GeV scenario. The displayed points 
comprise only a small fraction of the total number of points 
appearing in previous figures. 
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FIG. 41: A is plotted vs. k for the tan/3 = 1.7, msusY = 
500 GeV, A = -1000 GeV scenario. 

the couplings do not quite reach the non-perturbative 
value of coupling/ (47r) = 0.5 at the GUT scale. Rather 
coupling/ (47r) ~ 0.4 is a typical maximum value. In this 
sense they are the most attractive of the tan/3 =1.2 sce- 
narios. As for the blue points of this type in the tan/3 = 
1.7 case, the points in this special class at tan/3 = 1.2 typ- 
ically also have very modest A\, A^^ finetuning measure 
G, with G values between 10 and 30 being typical. One 
point of difference with tan /3 = 1.7 is that the tan /3 = 1.2 
special points all have B{hi — > aiai) > 0.75. 

For this smaller tan/3 = 1.2 value, B{t h^b) is 
larger 0.3) at the lowest mf^+ ~ 90 GeV mass than 
for tan/3 = 1.7. Nonetheless, the Tevatron is still unable 
to limit these scenarios since B{h^ T^Vr) < 0.2 (see 
Fig. US)) given the dominance of /i+ —> W^^*^b decays 
(Fig.SZl). 



III. CONCLUSIONS 

For low tan /3 values in the NMSSM, we have found 
many interesting new Higgs scenarios with a light CP- 
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FIG. 42: Ax is plotted vs. for the tan/3 = 1.7, msuSY = 
500 GeV, A = -1000 GeV scenario. 
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FIG. 43: cos^ 9a is plotted vs. rUa^ for the tan/3 = 1.2, 
msusY = 500 GeV, A = -1000 GeV scenario. 



odd scalar with mass below 10 GeV. For many of the ex- 
perimentally allowed parameter choices, the hi, /i"*" and 
/i2 are all sufficiently light as to be kinematically acces- 
sible at LEP and the Tevatron, but they decay into final 
states containing the light CP-odd scalar and therefore 
escaped detection. One particularly interesting parame- 
ter space region is that associated with the PQ symmetry 
limit (small |k| and small A^) in which the hi has mass 
near 90 GeV and the /i2 has mass near 100 GeV {i.e. ba- 
sically non-overlapping within experimental resolution) 
and g^2f^^B{hi — > bb) and gzzh2^^^^ ~^ ^) ^'^^ ^Vich as 
to explain the observed LEP excess throughout this re- 
gion. These points, such that both hi and /12 contribute 
to the LEP excess, are present for the tan/3 — 1.2 and 
tan/3 = 1.7 cases, but not for tan/3 = 2. 

Another important common feature of all these low- 
tan /3 scenarios that is also shared with the high-tan /3 
scenarios explored in earlier papers is that for any given 
TTiai there is always a lower limit on _B(T — > 701). This 
lower limit is above about 5 x 10~^ for < 7.5 GeV. 
We are hopeful that this is a level that can eventually 
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FIG. 44: is plotted vs. m^.j for the tan /3 = 1.2, msuSY = 
500 GeV, A = -1000 GeV scenario. 
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FIG. 45: is plotted vs. ruhi for the tan/3 = 1.2, mgusY = 
500 GeV, a'^ -1000 GeV scenario. 



be probed by BaBar and Belle. This lower limit arises 
because there is a lower limit on | cos Oa \ , and therefore on 
^aibb ~ cos0^ tan/3, below which B{hi — > aiOi) is not 
large enough for a light hi to have escaped LEP limits. 

We should further comment that all the scenarios of 
the present paper, as well as previous papers that fo- 
cused on higher-tan /3, are such that the scenarios that 
survive all experimental constraints are ones for which 
the contributions of the Higgs sector to both the pre- 
cision clectrowcak observables, AT and AS" (relative to 
the SM-Higgs contribution for m/ig„ = 110 GeV), and to 
the muon anomalous magnetic moment, (relative to 
the observed experimental discrepancy) are very small. 
In the case of AT and A^ the small extra AT can be 
understood as a natural result of either /13 and 02 decou- 
pling or of /i2 and 02 decoupling. 

Overall, the NMSSM provides a huge opportunity to 
have an "Ideal Higgs" boson scenario in which there 
is one or two light Higgs bosons (masses at or below 
100 GeV) that in combination have all the ZZ-Higgs cou- 
pling squared and therefore give values for the precision 



14 



tan/S=1.2, ;ti=150 GeV Msusy=500 GeV, At=-1000 GeV 
0.25 




0.00 



100 200 300 400 

m^. (GeV) 



500 



electrowcak observables 5* and T that are in excellent 
agreement with data. These Higgs bosons escape LEP 
limits because of unusual decays involving the light ai 
with iTLai < that is the common future of all these 
"Ideal" models. They also provide an excellent possibil- 
ity for describing the broad excess in the e+e" Zbb 
channel in the region mj^j e [90 GeV, 105 GeV] seen at 
LEP. We look forward to possibly discovering the ai in T 
decays at Babar or Belle, or direct detection of hi — > aiai 
at the LHC, if not the Tevatron. 
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